On the dependence of the reflection operator 
on boundary conditions for biharmonic 

functions 



Tatiana Savina a '* 



^Department of Mathematics, Condensed Matter & Surface Science Program, 
Nanoscale & Quantum Phenomena Institute, Ohio University, Athens, OH 45701, 

USA 



Abstract 

The biharmonic equation arises in areas of continuum mechanics including linear 
elasticity theory and the Stokes flows, as well as in a radar imaging problem. We 
discuss the reflection formulas for the biharmonic functions u(x, y) € K 2 subject to 
different boundary conditions on a real-analytic curve in the plane. The obtained 
formulas, generalizing the celebrated Schwarz symmetry principle for harmonic func- 
tions, have different structures. In particular, in the special case of the boundary, 
To := {y = 0}, reflections are point to point when the given on Tq conditions are 
u = d n u = 0, u = Au = or d n u = d n Au = 0, and point to a continuous set when 
u = d n Au = or d n u = Au = on Tq. 
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Introduction 



In this paper we study the dependence of the structure of the reflection oper- 
ator on different boundary conditions for the biharmonic functions, where a 
function u(x,y) of class C 4 (U) is said to be biharmonic if it is a solution to 
the equation A 2 u = [1]. Here U is a domain in M 2 , and A is the Laplacian. 
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The obtained reflection operator is generally an integro-differential operator, 
which reduces in the simplest case to the celebrated local (point to point) 
Schwarz symmetry principle for harmonic functions. 

In the case of harmonic functions, there are three basic types of boundary 
conditions: the Dirichlet, Neumann and Robin, and the Schwarz reflection 
principle can be stated as follows. 

Theorem 1.1 Let T C M 2 be a non-singular real analytic curve and P' G Y . 
Then, there exists a neighborhood U of P' and an anti-conformal mapping 
R : U U which is identity on Y , permutes the components U\,U2 of U \ Y 
and relative to which any harmonic function u(x,y) defined near Y and 

- vanishing on Y (the homogeneous Dirichlet condition) is odd HUtf . JF7[ /, 

u(x ,y ) = -u(R(x ,y )), (1.1) 

- subject to the Neumann condition onY, |^ = ; is even, 

u(x Q , y ) = u(R(x , y )), (1.2) 



- subject to the Robin condition on Y, a(x, y) 9u ^ y ^ +(3(x, y)u(x, y) = 0, (x, y) € 
T, can be continued by the following integro-differential operator |5]/, 

u(x , ?/o ) = u(R(x , y )) (1.3) 



R(x ,yo) 

J F{x,y,xo,y )u{u(x,y)) - u(x,y)u(F(x,y,x ,y )), 
r 

for any point (x ,y ) sufficiently close to Y. Here u(-) = -j^dx — -§^dy, and 
the integral is independent on the path joining an arbitrary point on Y with 
the point R(x ,y ). The function F is defined by the coefficients a(x,y) and 
(3(x,y), and the curve Y [5]. The mapping R mentioned above is given by 

R(xo,y ) = R(z )=Wo)i (1-4) 



where S(z) is the Schwarz function JB^- 

Note that if the point (x ,y ) G U\, then the "reflected" point R(x ,y ) G U 2 , 
and the mapping R depends only on the curve Y and is defined only near Y 
but may have conjugate-analytic continuation to a larger domain. 

The Schwarz reflection principle has been studied by several researchers (see 
PQ, 0, [T] [12] , [H]-[T8], [2l]-[27] and references therein). In particular, when 
T is a line, H. Poritsky proved that a biharmonic function u(x,y), a 
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solution to the biharmonic equation A£ w = 0, denned for y > and subject 
to conditions 

Ou 

u(x,0) = — (x,0) = (1.5) 



can be continued across the x-axis using the formula 

u{x , y ) = -u(R(x , y )) - 2y —(R(x Q , y )) - ylA x>y u{R{x Q , y )), (1.6) 

where R(x ,y ) = (x ,—y ) and A x<y = + J^. He also applied this for- 
mula to problems of planar elasticity, studying bending of plates, where u is 
a deflection of a thin plate, clamped along y = 0. An analogous formula has 
been obtained by R.J. Duffin [9j for three-dimensional case. Duffin also consid- 
ered spherical boundaries and applied his result to study viscous flows, among 
other things. J. Bramble [7] considered continuation of biharmonic functions 
across the circular arc with the clamped boundary conditions 

u= 9u =Q for ^ + 2 _ 2 qjn 

or 

He has shown that u can be continued using the formula 

u(x ,y ) = -u(R(xo,yo)) 

— ^2 — ( r o q^ u ( r ( x 0i Vo)) + ^( r o - a 2 )A x>y u{R{x , y ))), 



where r = \J x$ + y$ and a is the radius of the circle. J. Bramble also ap- 
plied it to the elastic medium problems in two- and three-dimensional space. 
Papers by F. John [H] and L. Nystedt [21] are devoted to further studies of 
reflection of solutions of linear partial differential equations with various linear 
conditions on a hyperplane. R. Farwig in [11J considered reflection principle 
for biharmonic functions subject to different boundary conditions on a hyper- 
plane. Following to H.A. Lorenz and R.J. Duffin [9], R. Farwig [11] applied 
his results to the Stokes system. 

The purpose of this paper is to derive and to compare reflection formulas 
for biharmonic functions across real analytic curves in R 2 subject to different 
boundary conditions occurring in a variational approach to A 2 . The motiva- 
tion of this study is as follows, most of the results for biharmonic functions 
mentioned above were obtained for the Dirichlet boundary conditions, see 
fll.5p and (jl.7p . However, unlike the Laplace equation, which is typically com- 
plemented with either the Dirichlet, Neumann or their linear combination, 
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the Robin conditions, the biharmonic equation has much larger set of possible 
boundary conditions as it is discussed in details in [13J . For this study we have 
chosen the conditions (i) u = d n u = 0, (ii) u = Au = 0, (iii) u = d n Au = 
(iv) d n u = Au = 0, (v) d n u = d n Au = and (vi) Am = d„Au = 0, appear- 
ing when one is considering a bilinear form associated with the biharmonic 
differential operator for u, v G C 4 (fl) [21] , [22] 

a(u,v) :=< A 2 u,v >= J vA 2 uduj. 

h 

Employing Green's formula to a(u,u), we obtain 

a{u } u) = J (u d n Au — d n u Au) ds + J(Au) 2 du. (1.8) 

an n 



The integral along the boundary dQ in (11.81) disappears under conditions (i), 
(ii), (v) and (vi), and therefore a(u,u) > 0. This immediately imply that 
the operator is strictly positive for (i) and (ii), and the corresponding bound- 
ary value problems in Q are well-posed. In the case (v), an extra condition, 

/ u duj = 0, is usually imposed. This yields the coercivity of the bilinear form 

n 

in suitable function spaces. 

Conditions (i), (ii) and (v) are commonly used in applied problems for the bi- 
harmonic equation as well as for other equations with the biharmonic operator 
in the principal elliptic part. To give some examples, we remark that condition 
(i), the Dirichlet condition, corresponds to the clamped plate model [7], [9], 
[25] . while condition (ii), the Navier boundary condition, corresponds to the 
hinged plate, when the contribution of curvature is neglected [B], [13]. In the 
later case the plate is ideally hinged along all of its edges so that it is free 
to rotate and does not experience any torque or bending moments about the 
edges. Both the Dirichlet and Navier conditions are also used to model elec- 
trostatic actuation |20j alone with the governing non-linear nonlocal elliptic 
equation with the bi-Lanlacian in the principal part. 

Condition (v) for biharmonic equation was considered in |28j, moreover, this 
condition is often used for the famous Cahn-Hilliard equation, u t = — A(Ati + 
u — u 3 ), which is a semi-linear parabolic equation describing (among other 
processes) spinodal decomposition [23], [28]. Pattern formation resulting this 
phase transition (spinodal decomposition) has been observed in alloys as well 
as polymer solutions and glasses. In this case the physical meaning of the 
second condition in (v), d n Au = 0, is a no-flux condition (none of the mixture 
can pass through the walls of the container). The first condition in (v), d n u = 
0, is the most natural way to ensure that the total free energy of the mixture 
decreases in time (this is a requirement from thermodynamics - the variational 
condition). 
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We have to comment on condition (vi), for which we will not be able to derive 
a new reflection formula. Indeed, consider biharmonic function in U subject to 
conditions (vi) on T. Denote v(x, y) := Au(x, y), then function v is a harmonic 
function in U subject to conditions v = d n v = on T. Thus, v = in U, and 
therefore u(x,y) is a harmonic function in U. Since any harmonic function is 
a solution to the problem in question, it means that the space of solution does 
not have finite dimension. Note also that as it was shown in [13], condition 
(vi) does not satisfy the Complementing Condition, introduced by S. Agmon, 
A. Doughs, and L. Nirenberg [2J. This condition is necessary for obtaining 
estimates up to the boundary for solutions of boundary value problems to the 
elliptic equations and, therefore, is crucial for the existence and uniqueness 
results [2J. 

In the case of conditions (iii) and (iv) only one of the term in the boundary in- 
tegral (jl.8p disappears, so our interest to this conditions is mostly theoretical, 
however, one of the problems arising in radar imaging may be reformulated 
as a boundary value problem for the biharmonic functions with condition 
d n u = u — Au = 13], which coincide with (i) if (3 = and with (iv) if 
(3 = — oo. We remark, that other linear combinations, the Robin-type condi- 
tions, involving (i)-(vi) are used in applications as well, but the corresponding 
reflections will be discussed elsewhere. 

The reflection formula for the case of condition (i) was obtained in [TJ. The 
only known in the literature result for the conditions (ii) and (v), is in [TT] . 
where the boundary conditions are given on x-axis, and the odd point-to-point 
reflection holds for (ii) and even for (v). The author is not aware of any results 
for the conditions (iii) and (iv). Our aim is to derive reflection formulas for the 
conditions (ii)-(v), that is, the formulas expressing the value of a biharmonic 
function u(x,y) at an arbitrary point (xo,2/o) Ui in terms of its values at 
points in U2, when the boundary conditions are given on a real-analytic curve, 
and to study the properties of these mappings. 

We remark that the structure of the reflection formulas attracts attention 
of many researchers, interested, in particular, to answer the question: when 
the reflection is point to point [15], [10], point to final set [18J or point to 
continuous set [5]. 

The structure of the paper is as follows. In the section 2 we discuss the so-called 
reflected fundamental solution for each case, whose properties depend on the 
boundary condition and, therefore, determine the structure of the reflection 
operator. In the section 3 we formulate and prove the main theorems. 
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2 The reflected fundamental solutions 



This section is devoted to one of the key steps in deriving the reflection for- 
mula, that is, to the construction of the reflected fundamental solutions for 
each case of the boundary conditions. 

Let r C M 2 be a non-singular real analytic curve and point P' G V . Consider 
the biharmonic differential operator in a neighborhood U C M 2 of the point 
P' . Let U\ and U2 be the components of U \ T, and P(xq, yo) be a point in U\ 
sufficiently close to V. The fundamental solution can be written in the form 

G = -t^-((> - xo) 2 + (y~ Vo) 2 ) ln((ac - x ) 2 + (y - y ) 2 ) + g{x, y, x , y ), 

1D7T 

where g is a regular biharmonic function. It is obvious that G is a real-analytic 
function in M. 2 except for the point P(xq, yo) G U\. 

Let us complexify the problem, that is, consider a complex domain W in the 
space C 2 to which the function / defining the curve V := {f(x, y) = 0} can be 
continued analytically such that W n M 2 = U. Using the change of variables 
z = x + iy, w = x — iy, the equation of the complexified curve r c can be 
rewritten in the form 

(z + w z - w\ _ 

and if grad f(x,y) 7^ on T, can be also rewritten in terms of the Schwarz 
function and its inverse, w = S(z) and z = S(w) [8]. In the characteristic 
variables G(z, w, z , w ) can be rewritten as 

G = -—[-((* - zq)(w - w ) ln[(z - z ){w - w )} + g{z,w, z ,w ). (2.9) 

l07T 



It is obvious that the continuation of G to the complex space (12.91) has loga- 
rithmic singularities on the complex characteristics passing through this point, 
i.e., on K P := {{x - x ) 2 + (y - y ) 2 = 0}. 

Note, that the specific choice of the regular part g of the fundamental solu- 
tion G does not affect our final result. Thus, for convenience we choose the 
fundamental solution in the form 



G(z, w, z , w ) = --— w, z , w ) + G 2 {z, w, z , w )), 

lD7T 

G x = (z- z )(w - w )(Hz - z ) - 1), (2.10) 
G 2 = (z - z )(w - w )(\n(w - w ) - 1). 



Our goal is to construct (multiple-valued) functions j = 1,5, which are 
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biharmonic functions satisfying on Tq one of the pair of conditions (i) G — G 
and dnG = d n G, (ii) G — G and AG = AG, (iii) d n G = d n G and AG = AG, 
(iv) G = G and <9 n AG = <9„AG, (v) d n G = d n G and d n AG = <9„AG, and 
having singularities only on the characteristic lines intersecting the real space 
at point Q = R(P) (see formula (11.41) ) in the domain {7 2 and intersecting Tq 
at Kp D Tc- These functions are called the reflected fundamental solutions. 
According to (I2.10p it is convenient to seek these functions in the form 

G i:,) {z,w,z ,w ) = -—-(G^\z,w,z ,w ) + G [ 2 3 \z,w, z ,w )). (2.11) 

l07T 



It is easy to check that for the case (i) functions Gi and G^ are 

Gi = (z - z )(w - w )(ln(S(w) - z ) -l) + (z- S(w))(w - w ), 
G« = (z- z )(w - w Q )(ln(S(z) - w ) - l)+(w - S(z))(z - z ). 

Similar functions (up to the regular part) were used in [TJ, where function g 
in the expression for G (12. 9p was different from the chosen in (I2.10p . 

To obtain G± and G% for the cases (ii) - (v), we use the asymptotic expan- 
sions introduced by Ludwig [151 . seeking these functions in the form 



oo 

G? = J2 b k\ z , w , w o)fk{S{w) - zq), (2.12) 



Gf = Y,af{z,w,z G )f k {S{z)-w Q ), (2.13) 
fc=i 



where 



/A(0 = |j(lne-Cfc). k = 0A, 



fk(0 = (-iy k ~ l (-k - l)!£ fe , k < -1, Go = 0, (2.14) 



k x 



Gk = Y.~ v fc = l,2,.... 



i=i 



This form ensures the location of singularities on the desired characteristics. 
Substituting ( 12. 12ft and ( 12.131) into the biharmonic equation, we obtain the 
following differential equations for the coefficients and 



d 2 a^ ] d 2 b[ j) 



dw 2 dz^ 



1,2,.... (2.15) 



7 



To find coefficients af* and for each case we use the conditions on T c . 
Case (ii), = G and AG (2) = AG on T c , 



a® = z-z , ^S' = l onr c , 

ow 

rtn® f) 2 n® 

OW owoz 



hi 2 ) db ? ] q> i r 

= w — u> , — ^ — 5 =1 on T c , 



Thus, 



dz 

K "U, ^ & - dwdz , k-2,6,.... 



a? : :„-^^. «« (.,)-^^/;-. (2.16) 
6^= W - Wo + (g ~l (w)) , ff> = (-l)^ H^ Dr, (2.17) 



where D 2 = 4 2 = etc, and D w = 

Note that for the special case when T is a straight line y = 0, we have 

GS 2) = (z- w )(w - 3b)(ln(u> - zo) - 1), (2.18) 
G^ 2) = (w- z )(z - w )(\n(z - w ) - 1). (2.19) 

Case (m), d n &V = d n G and AG (3) = AG on T c , 

a® = -(z-z ), ^t S ' = 1 onr € ? 

(3).,, da® da® 0a® d 2 a® 

a 2 O — 1 + D , O — rv rv 5 

oz OW OW OWOZ 

a (3) s , = da® x da® l g/ da® gl = d 2 a® 1 fc > g 

fc dz dw dw dwdz ' 



ai,(3) _ 
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dw dz dz dwdz 1 

ql(3) ™(3) ™(3) ^2l(3) 

^(S)^/ = dt) k~l + 0°k-l gl dt) k g, = g °k-l k>3 

k dw dz dz dwdz ' — ' 

.f = -(* - ,0) + bf = -(„ - W „) + <^|M, (2.20) 



(3) 



(-1)*(|4*- 2 - ^ s} Z)) Dl-% k>2, (2.21) 
6W = ( _ 1) *(^2_(iz|H)^j r i) > fc > 2) (2.22) 



where D° = D° = 1. 

For the special case when T is a straight line y = 0, we have 

(5^ = — (2w — z — wq)(w — z )(ln(w — zq) — 1) 

+2(w - z ) 2 (\n(w - z ) -3/2), (2.23) 

= — (2-2 — U> — 2o)(z — U>o)(rn(z — Wo) — 1) 

+2(2- w ) 2 (ln(^- w ) -3/2). (2.24) 
Case fwj, G (4) = G and <9 n A<5 (4) = <9 n AG on r c , 



Q (4) 

a[ 4) =z-z , ^S' = -l onr c , 

a„(4) ~2 n (4) a (4) 

dw dwdz dw 

a (4) =Q? ggfc (5") 2 = 2 fc ~ 1 5 / aa ^ s" ^ a fc-2 ; A , >3; 
fc dw dwdz dw dwdz 2 ' ~~ 



&i = w - wo, = -1 011 r c, 

f)/,(4) ^2l(4) or (4) 

aw dz 
owaz dzdw 2 
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of = <^>>, af = (-If { -^p£»b*-\ (2.25) 

b W =w _ WQ _ (izjWl , f,W = * > 2.(2.26) 



For the special case when T is a straight line y = 0, we obtain 

G{ 4) = (2w-z- w )(w - z )(ki(w - z ) - 1), (2.27) 
G { 2 ] = (2z-w- z )(z - w )(\n(z - w ) - 1). (2.28) 

Case (v), d n G^ = d n G and d n AG^ = d n AG on T c , 



<9a (5) 

4 5) = -(z - z ), ~dtj~ S ' = ~ l ° n Fc ' 

f) n (5) a_(5) o (5) «2 (5) a (5) 

(5) c /_-, ggl , aa l q/ /o/\2_ q ' Q l q/ o// 

ot ot cw aw 

aQ fc /gA2_ 2 S' Ua k-\ cj/l ° a k-2 ^>3 

<9w dwcfe <9u> dwdz 2 ' 



<% (5) ~ 

6i 5) = -(w - w ), ~{t~ S ' = _1 ° n r<C ' 

owoz oz 

, (5) 5/ _ dfrfc^i , dfrj-i q/ 

dz dwdz dz dzdw 2) 

af = _( 2 - Zo) - <HLZ*M>, (2 . 29 ) 



-D'l^fl!--^ 1 ^'), *>2, (2.30) 



a (5) 



l P> = -(„-^)-(iz|W), (2.31) 
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(2.32) 



For the special case when T is a straight line y = 0, this imply 

Gj 5) = -{z - w )(w - z )(Hw - z ) - 1), (2.33) 
= ~(w - z )(z - w )(\n(z - w ) - 1). (2.34) 

Thus, we have constructed the reflected fundamental solutions for each case 
of the boundary conditions (ii)- (v) as formal series. Convergence of the series 
follows from convergence of the multipliers of the logarithms in (I2.12p and 

Vl»=p}t\z,w, Wo f S{w) - Z <> )k , (2.35) 

V? = t-'i\^) (S(Z) Z W0)k - (2-36) 

k=l 



Coefficients in ( I2.36P do not depend on w$, therefore, this expression can 
be interpreted as the Taylor series of V 2 U) function of — w at the point 
-S(z). Function V^ j) can be described as the unique solution to the Cauchy- 
Goursat problem for biharmonic functions with holomorphic data on Tc and 
the characteristic line S(z) = wq. The data is prescribed by the boundary 
conditions for each j = 2, 5. For example, for j = 2 the corresponding problem 
is 



in W 

on r c , 

on T c , (2.37) 
on S(z) — wo = 0, 

on S(z) — Wq = 0. 

Note, that S(z) is an analytic function in the neighborhood of Tq and its 
derivative does not vanish on T^. Existence and uniqueness of holomorphic 
solutions to the Cauchy and Goursat problems for holomorphic partial differ- 
ential equations with holomorphic data are discussed in [16]. 

The reflected fundamental solutions as well as functions Vi and Vj; are used 
in the next section for deriving the corresponding reflection formulas. 



dz 



V, 



(2) 



A 2 Vf ) = 0, 
(z - z )(w - w ) 
AV 2 (2) = 1 




(2) 



(z - z )S' + (w- S(z)) 
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3 The main result 



First we state the reflection formulas for biharmonic functions given in the 
upper half plane describing the analytic continuation across the the x-axis. 

Theorem 3.1 Let U cR 2 be a sufficiently small domain divided by a straight 
line r C U := {y = 0} into two components U\ C and U 2 C M?_ ; 
(xo,yo) G U\ and (xo, —yo) £ LJ%. Then any biharmonic junction u(x,y) in the 
domain U\ subject to one of the conditions (i)-(v) on Tq can be continued to 
the domain U2, using the following formulas: 

(i) If u = d n u = on T , then 

du 

u(x , y ) = -u(x , -y ) - 2y — (x , -y ) - y^A Xty u(xo, -y ), (3.38) 

(ii) if u = Au = on Tq, then 

u(x , y ) = -u(x , -y ), (3.39) 

(Hi) if u = d n Au = on To, then 

-yo 

u(x , y ) = -u(x , -y ) - y J Au(x , y) dy, (3.40) 



where the integral is computed along the segment parallel to y-axis, 

(iv) if d n u = Au = on T , then 

-yo 

u(x ,y ) = u(x ,-y ) - J yAu(x ,y)dy, (3.41) 



(v) if d n u = d n Au = on To, then 

u(x , yo) = u(x , -y ). (3.42) 

Remark 3.1 // a biharmonic function u(x,y) is also a harmonic function, 
then Ii3. 39\) and Ij3.40\ ) coincide with the odd continuation U.l\) . while formulas 
an d &3-4ty) with the even continuation hi. 2^1 for harmonic functions. 



Theorem 3.2 Let U C M 2 be a sufficiently small domain divided by a non- 
singular real analytic curve T into two parts U\ and XJi- Let also P(xq, yo) be a 
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point in JJ\ having its reflected point Q(R(xo,yo)) in U2 (see Then any 

biharmonic function u(x,y) in the domain U subject to one of the conditions 
(i)-(v) on the curve V C U can be continued across T, using the following 
reflection formulas: 



(i) if u = d n u = on T, then 



u(P) = -u{Q) - [x 



2 > dx 

2i hy' 



/ S(x + iy ) - S(x -iy )\du 
" ( y ° + m ) ThiW 



-\i x l + Vo~ s i x o + m){xo + iy ) (3.43) 
-S(x Q - iy ){x - iy ) + S(x + iy )S(x - iy ))A x>y u(Q), 

(ii) if u = Au = on T, then 

u{P) = -u(Q) + K 2 , (3.44) 

(Hi) if u = d„Au = on T, then 

u(x , y ) = -u(x , -yd) + K3, (3.45) 

(iv) if d n u = Au = on T, then 

u(x Q , y ) = u(x Q , -y ) + K 4 , (3.46) 

(v) if d n u = d n Au = onT, then 

u(x , y ) = u(x , -y ) + K5, (3.47) 

where 



1 ? („uydAu , dV& >„,*du dAV^' 



^ = T- / Vi3) ^T- - Au ^T~ + AV(J) 7T - ^ n dx (3.46 

8i j \ ay ay ay ay / 



r 



- v a Am ^T~ + Ay V ~ ^ o rf y, 

y ax ax ox ox J 

the integral is computed along an arbitrary path joining the curve V with the 
reflected point Q, V® = V} j) - V 2 U) , j = 275. 
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It is obvious that the Theorem 3.1 is a special case of the Theorem 3.2., thus 
we will prove only the later for the boundary conditions (ii)-(v) (for the case 
(i) see HI). 

For simplicity, we assume that T is an algebraic curve. Under this assumption, 
the Schwarz function and its inverse are analytic in the whole plane C except 
for finitely many algebraic singularities. 

Proof of the Theorem 3.2. The main step of the proof is already done by con- 
structing the reflected fundamental solution for each case of boundary condi- 
tions (see section 2). The rest of the proof is based on the contour deformation 
in Green's formula [12] and is similar to pQ. 

The Green's formula, expressing the value of biharmonic function at a point 
P via the values of this function on a contour 7 C U\ surrounding the point 
P, is 




<P) = I ( G ^ ~ A -f7 + ^ - u^pj dx (3.49) 

. dG A Bu dAG\ , 
- Au— + AG— - u — dy, 

ox ox ox J 

where G = G(x,y,x ,y ) is an arbitrary fundamental solution of the bi- 
Laplacian. The most suitable one for what follows is (I2.10p . 

Since the integrand in (I3.49P is a closed form, the value of the integral does 
not change while we deform the contour 7 homotopically. Thus, the goal is 
to deform the contour 7 from the domain U\ to the domain U2 by deforming 
it first to the complexified curve Tq. This part of the deformation is possi- 
ble if the point P lies so close to the curve T that there exists a connected 
domain Q C Tq such that Q contains both points of intersections of the 
characteristic lines passing through the point P and Q can be univalently 
projected onto a plane domain (for details, see [2Q\). Thus, we can replace 
the contour 7 in ( I3.49|) with the contour 7/ C Q, which is homotopic to 7 in 
C 2 \ {(x - * ) 2 + (y - y ) 2 = 0} =: C 2 \ K P . 

Note that due to homogeneous boundary conditions (ii) - (v) a half of the 
terms (different for each case) in the integrand of (I3.49p . while integrating 
along 7/, vanishes. For example, formula (I3.49P in the case (ii) can be rewritten 
in the form 

<P)=l{ C^ + M%)*-(G^ + **)«. (3.50) 
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To deform the contour 7/ from Tc to the real domain U2 we replace the fun- 
damental solution with the corresponding reflected fundamental solution G^\ 
j = 1,5 (see section 2). Since functions G^> have singularities only on the 
characteristic lines intersecting the real space at point Q = R{P) in the do- 
main U2 and intersecting Tc at Kp D Tq, we are able to deform contour 7/ 
from the complexified curve Tc to the real domain U2 without changing the 
value of the integral [26]. As a result, we obtain 



« P = / G w — Aw—- + AG W - - it — — dx 

dy dy ay ay J 



7 



~ G "a Aw^— + AG W_- u _ dy, (3.51) 

y ax ax ax ax / 

where 7 C C/2 is a contour surrounding the point Q and having endpoints on 
the curve V. Formula (I3.5ip in the characteristic variables has the form, 



}/ , , r , (; , . ■■ | d 2 G& du u d 3 G^ d 2 u 9&^ \ dz 
.•1- ■».•)... Q^^u; <9z dz 2 dw dzdw dz ) 




+ d 2 G^ du _ u d z G^ _ d 2 u dG^ \ ^ g 
dzdw 2 dzdw dw dzdw 2 dzdw dw J 

If we substitute (12.111) . (I2.12p and (12 . 13[) into (13. 52|) and move one endpoint 
of the contour 7 along the curve T to the other endpoint, the integral terms 
containing products of the function u and regular part of the function G^' 
and their derivatives vanish, while the integral terms containing logarithms 
can be combined as follows, 



— fln(S(w) z )Uv U) I d2V ^ du ° 2u dV ^ )dz 

An J W I 1 dz 2 dw dzdw dz dzdw dz 

7 

(V U) + _ U FVP _ d 2 u , 

dzdw 2 dzdw dw dzdw 2 dzdw dw > ' 



— / hs(z) - w )\(v U) d ' u + ^2_^ _ u &Yll _ d ' u dv *\ dz 

Att I \ 2 dz 2 dw dzdw dz dz 2 dw dzdw dz 

7 

-(VP— + - — ^)dw\, (3.54) 

dzdw 2 dzdw dw dzdw dw > 
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where 7 is a loop surrounding the point Q and having endpoints on the curve 

r. 

The logarithm ln(5(w) — zq) obtains the increment 2ni along the loop, while 
the logarithm ln(S(z) — Wq) gets (— 2iri). Thus, compressing 7 to a segment 
joining Q to T, we find that the integral (13.531) can be rewritten as 



1 f[ (v g) d 3 u d 2 V^ du d 3 V 2 (j) d 2 u dV& . 

2 J dz 2 dw dzdw dz dz 2 dw dzdw dz 



(j) d 3 u | d 2 V^ du u d 3 vl j) Q2u dy{j) ) dw \ (355) 
dzdw 2 dzdw dw dzdw 2 dzdw dw ' ' 

where V® = v} j) - V 2 (j) . 

The rest of nonzero terms in (13.521) are terms involving derivatives of the log- 
arithms. Some of these integrals also vanish due to the properties of functions 
Vi and V 2 \ resulting in 



r, da? (S'(z)) 2 u dz db? ((S'(w)) 2 u ^ 

47riv dw S(z)-w Z+ dz S(w)-z W ) 

7 

1 dn (j) 

= -2 um ~dt {Q)s ' {Q) + -^(Q) S '(Q))- ( 3 - 56 ) 

Combining (" 13 . 55[) and (I3.56P we finally obtain, 
1 dn {j) dh^ 

u(P) = -- U (Q)(^-(Q)S'(Q) + ^-(Q)S'(Q)) + K r (3.57) 



Here the expression in the parentheses with the appropriate choice of the 
coefficients and b± (see section 2) is equal to either 1 or -1. Thus, formula 
(I3.57P in variables x,y is equivalent to (I3.44p - (I3.48p . 



Remark 3.2 Formula (3.51) gives continuation of a biharmonic function 
from the domain U\ C IR 2 to the domain U2 C IR 2 as a multiple-valued function 
whose singularities coincide with the singularities of the functions S or S . 

Remark 3.3 In the special case when T is a straight line and boundary con- 
ditions (ii) or (v) are applied, = V%\ and, therefore, Kj = 0, j = 2,5. 



Example 3.4 As an example of formula flff.JTp consider a biharmonic func 



tion u(x,y) subject to the Navier conditions (ii), u = Au = 0, on a unit circle 
centered at the origin, x 2 + y 2 = 1. In this case series 112. 35\) and \2.36) with 
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coefficients A2.16]) , can be summed, and = (z — z )(l/z — w ) + 



[1/z — w)(1/wq — z). The reflection formula then has the form 



u(r , 6» 



1 r 2 — 1 

r 4r 



(~<9 r u(r, 6q) + d 2 rr u(r, j) dr, (3.58) 



where the integral is computed along the straight line 6 = 9q. 

Acknowledgments. I would like to thank Professor Reinhard Farwig for 
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